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Capítulo 1: Operadores em Química Quântica 


1.1 Introdução 


A Mecânica Quântica foi construída na década de 20 para lidar com o átomo, os 
físicos consideravam-na como um instrumento provisório destinado ao fracasso fora dos 
domínios do átomo, entretanto a mecânica quântica prosperou e ultrapassou os mais 
fantásticos sonhos de seus inventores. Atualmente, devido a aceitação dos Postulados 
Quánticos e suas extensivas aplicações em Centros Científicos, vem proporcionando maior 
divulgação de experimento clássicos da Mecânica Quântica, bem como seus inventores. A 
internet, onde é um meio de informação de mais alto padrão, nos favorece diversas fontes 
de trabalhos diretos ou correlacionados à Mecânica Quântica, sendo estes: revistas, 
fascículos e até mesmo jornais científicos. 


1.2 Conceitos Preliminares 


Para explicar a mecânica dos átomos е dos sistemas nucleares, foram desenvolvidas, 
a partir de 1900, diversas teorias. Estas teorias foram se enquadrando a uma denominação 
comum: “Mecánica Quántica" (MQ). 

Assim, MQ é um conjunto de teoria dos sistemas atômicos e nucleares. Este conjunto 
não é unitário, mas sim uma sucessão de diversas teorias, umas complementares das outras. 
Estas teorias surgiram e desenvolveram-se da Física Clássica, particularmente da mecânica 
newtoniana e da teoria eletromagnética de Maxwell. Demonstrando-se que os conceitos 
clássico e eletromagnético eram suficientes para explicação do mundo macroscópico, mas, 
incapazes para uma explicação adequada e coerente do mundo microscópico, i.e., atômico- 
nuclear e também molecular, nasceram à chamada “Antiga teoria dos quanta” ou Mecânica 
Quântica Antiga. Esta mecânica compreende as teorias de Planck (1900), Einstein (1905), 
Bohr (1913) e De Broglie (1924). 

A antiga MQ embora explicasse muitos fenômenos até então incompreendidos, 
falhava em sua base lógica para outros fenômenos. E assim Heisenberg (1925), 
desenvolveu a “Mecânica das Matrizes”, que não implica em nenhum modelo atômico, pois 
é um edifício puramente matemático. Quase simultaneamente, Schrôndiger (1926) 
baseando-se nas idéias de De Broglie, desenvolveu a “Mecânica ondulatória”. Esta 
mecânica tem um aspecto mais físico, embora algo impreciso. Esta imprecisão foi 
desaparecendo com a aplicação das teorias de Pauli (1924) e Dirac (1926). As teorias de 
Heisenberg até Dirac, são conhecidas como “Mecânica Quântica Moderna”. 
Esquematicamente, temos: 


teoria dos quantas de Planck 
MQA | teoria atômica de Bohr 


dualismo de De Broglie 
MQ 9 


MQM | mecânica das matrizes de Hiesenberg 


mecânica ondulatória de Schôndiger 


teoria do Spin de Pauli 


teoria da energia positiva e negativa de 
Dirac 


Juntamente com estes físicos, representantes da МО, figura outro grande número de 
cientista que, experimentalmente, deram confirmação a cada um destas teorias, assim como 
aplicaram teoricamente a muitos fenômenos Físicos-Químicos. 

Primeiramente estaremos considerando a visão matemática, que a mecânica 
quântica possui com uma linguagem acessível para estudantes de ciências e assim uma 
breve discussão da MQA (Planck de Broglie) e MOM (Pauli Dirac). tratando 
historicamente os fatos ocorridos e descrevendo alguns passos tomados pelos pesquisadores 
da MQ, explanando o formalismo matemático. 


1.3 O Formalismo Matemático da Mecánica Quântica 


Na mecánica quântica as formas de se estabelecer relações numéricas е 
representação física de eventos são de modelos não convencionais como o de costumeiro da 
física e química. As vezes tendo tratamento até filosófico de problemas. 


1.4 Operadores 


Em muitos processos matemáticos um elemento tem que se submeter a um conjunto 
de operações quase sempre complicadas, mas que se repetem muitas vezes, a fim de aplicar 
a outros elementos diferentes para correlacioná-los. O conjunto de operações e 
transformações métricas são representadas por símbolos (À, V?, П, V) que se denominam 
operadores de transformação. Este conjunto de operações de transformação pode ser 
representado também por uma matriz, a qual tem propriedades intrínsecas da teoria de 
matrizes. Veja abaixo a representação de operadores e matrizes respectivamente. 


f = Ah [1] 


А= а dy 4; 2] 


Estes operadores são denotados usando letras com acento circunflexo, como por 
exemplo Ô. Existe em princípio, um operador para cada informação que se deseja sobre o 
sistema, sendo os mais importantes os operadores que fornecem a energia, chamados de 
operadores Hamiltonianos que são denotados pelo símbolo H. Os operadores, por serem 
ordens ou instruções matemáticas, devem ser aplicadas pela direita nas funções desejadas 


por exemplo a função de onda Y ou usando a chamada representação de Dirac ЕШ Se Ô = 


d/dx, ou seja a ordem matemática é “derive com relação a x”, e Y = f(x),podemos ver 
claramente que Ó | JP) (df(x)/dx) não é o mesmo que |Ҹ)Ӧ , ou f(x)d/dx (4). Como em 


alguns casos é necessário trabalhar com combinações lineares de funções de onda, os 
operadores também devem ser lineares, 1.е.: 


0(%,) + |w,)) - О|Ф,)+ Ó| v) [3] 


Quando aplicamos о operador na função de onda pode acontecer dois eventos. No 
caso mais simples a operação sobre V fornece a própria função V multiplicada por uma 
constante a, ou seja: 


Ó|W) = a|W) [4] 


Neste caso a função de onda Y é dita ser autofunção do operador Ô com autovalor 
a. O autovalor obtido corresponde à propriedade que se deseja medir, de forma que, quando 
se fizer uma medição desta propriedade em um conjunto de sistemas idênticos, sempre será 
obtido o mesmo valor a. Mas é óbvio que as medições realizadas usando os operadores 
devem fornecer valores reais e não números complexos. Apresentando esta condição o 
operador é dito Hermitiano, o que matematicamente corresponde a expressão: 


(“ДО ғ)-( |Ó v) [5] 


А equação [5] garante que se Ô| +) = а +) e que дз” = aw) ¿entšoa =a ‚о 


que só acontece se а for um número real, ou seja, se não contiver o número i na sua 
expressão. 

O segundo tipo de resultado que pode ser esperado quando se aplica um operador 
em uma função de onda é que a função de onda seja alterada pelo operador, como mostrado 
na equação 6: 


Ô| Y) =b) Y) [6] 


Agora a função não é autofunção do operador, e o valor b corresponde a uma única 
medição da probabilidade desejada em um só sistema dentro de um universo de sistemas 
idênticos. Neste caso , cada vez que se fizer a medição experimental da propriedade seria 
obtida uma resposta diferente, sendo portanto necessário calcular o valor médio das 
medições, também chamado de valor esperado, <b>, para ter um resultado confiável da 
medição. O valor médio é calculado exatamente como se calculam os valores médios em 
estatística: 


(v Ó|v) 
V) Gy а 


O principal problema agora é definir matematicamente os operadores е as funções 
de onda. Os operadores são obtidos através dos comutadores dos operadores. 
Se tivermos uma função de onda descrita por 


2лі 
E (ap, Et) 


Y = Ae ^ [8] 


que pode ser aplicada а descrever qualquer partícula de massa т е velocidade У, que 
apresente comportamento ondulatório, sendo portanto uma equação geral e não só referente 
ao fóton. Agora utilizando-se dum operador para derivá-la em relação à x, e considerando 
que a energia E e o momento P, da partícula associada à onda são constantes, temos: 


[9] 
rearranjando temos: 


к y [10] 


analisando a equação [10] vemos que a função de onda Y é uma autofunção do operador 
ih O р 

ш уал» com autovalor p, , ou seja, que о operador extrai o valor de momento da função 
л Ox 


ih O € š 
de onda — DUE Isto significa que, para obter operadores podemos escrever a expressão 
л Ox 


clássica da quantidade que nos interessa, e substituir os termos referentes a componente em 


ih O ES 
x do momento, р,, por — 225 Por analogia, para as componentes em y е z a substituição 
л OX 


ih O Š ih 
2л Оу 2л 02 
tempo (t) não é necessário fazer nenhuma modificação para converter a expressão em um 
operador. Esta metodologia pode ser expressa em forma de tabela. Veja Tabela 1 


deve ser por — respectivamente. Já onde existem coordenadas (х,у,2) ou 


Tabela 1: substituições para a construção de operadores. 


Termo na Termo no Símbolo do 
expressão clássica operador operador 
X X $ 
Y Y y 
Z Z 2 
Р, ih O P. 
2л Ox 
Р, — ih 0 ГА y 
27 Oy 
P, ih 0 P. 
ETT 
T t г 


Outros exemplos que expressa estes operadores pode serem vistos nas Tabelas 2 e 3 


Tabela 2: Exemplos do operador Hamiltoniano para о movimento de uma 
partícula de massa m em diferentes campos de força definidos pela função 
(operador) potencial V. 


Operador 
r р 2 2 
(а) Partícula livre .. v0 0. f. M 
x 2m dx? 
b) barreira de ^ 
(b) V=V, H= Ls (х<0; x>0) 
potencial PEN 2m d 
» 2 2 
0 a x Н- -— а +V, (0<х<а) 
А 2 2 
© " Oscilador A V=(1/2)K2 p. h d lge 
armónico 2m dx? 2 


Tabela 3: Exemplos do operador para átomos e moléculas. 


Operador H 
(a) Atomos de Bor n? - n Ze? 
um elétron 77 am ЖЕШ 


(b) Átomos de 


Ze 
muitos elétrons уе D : XXL 


i j»i "ij 
(c) moléculas po "mese Z, LL 
ча. "т E res : хуу : 
гі i j>i 17 а В>а Tap 


1.5 Comutadores 


Em 1926, Born, Heisenberg e Jordan apresentaram na Zeitschrift für Physik a 
mecânica quântica, desenvolvida por Born em 1924, e Heisenberg, em 1925, sob o ponto de 
vista matricial. Nesse trabalho, eles apresentaram as relações de comutações para o 


momento angular L para um sistema de muitas partículas: 


h 
БЕДЕ a É 


h 
IL.L.]- 7L. [11] 


h 
МШБ b. 


Hoje em alguns cursos de MQ, podemos ver que esses comutadores podem ser visto 
da forma apresentada acima. 


Mas qual será a idéia de comutação 2, comutar significa “Realizar diversas 
funções”. Imaginemos uma situação na qual temos dois operadores Hermitianos Ó e É que 
comutam entre si, ou seja ÔÊ — ЕО = 0 (o comutador também é representado como [Ô,É] 
= 0). Suponhamos ainda que existe um conjunto de funções de onda V; que são autofunções 
do operador Ó, com autovalores diferentes entre ai, i.e., sem degenerescéncia”. 


O|Y,)=aY,) 


neste caso qual seria a relação do operador E com este conjunto de funções?. Se operarmos 
pela esquerda com É teremos que 


ÉÓ|v.) = E(a PP) = а(8|%,)) [12] 
e pela condição de comutação: 
ЕО = ОЁ > ÉÓ|w) = О(Ё|\Р,)) [13] 
de (12) е (13) temos que 
О(Ё|\Р,) = а,(Ё|\Р,)) [14] 


A equação [14] nos mostra explicitamente que É | ,) é uma autofunção do operador 


О, com o autovalor а. Como V; é a única função com autovalor a; , devido “a nào 
degenerescência dos sistema, a única forma da equação [14] ser verdadeira é quanto a 


função Е|ч) seja um múltiplo de V , i. e., 
Ё|\Р,) = aw) [15] 


Em outras palavras, а função |е) também é autofunção do operador É. Este 


resultado na forma de teoria seria: “Se dois operadores Hermitianos О e É comutam, então 
existe um conjunto de funções que são autofunções dos dois operadores”. Este teorema é 
extremamente importante pois mostra que as únicas propriedades de um sistema que 
podemos medir simultaneamente e com precisão são aquelas para as quais a função de onda 
representando o estado do sistema é a autofunção dos seus respectivos operadores, ou seja, 
para medir precisa e simultaneamente duas ou mais propriedades de um sistema, seus 


* Degenerescência: significa os mesmos autovalores de energia. 


respectivos operadores devem comutar. De fato, este resultado é uma forma de expressão 
do Princípio de Incerteza de Heisenberg, pois indica que só é possível conhecer 
simultaneamente com precisão as propriedades de um sistema cujos operadores comutam. 
Claramente, de acordo com a formulação mais conhecida deste Princípio, os operadores 
referentes à posição e ao momento de um elétron em um átomo não comutam. 


1.6 Observáveis 


As propriedades (E, P, Х;) que formam um observável têm duas características 
importantes: 

1) Podem excluir-se mutuamente, ou 2) Constituem um conjunto completo, 1.е., incluem 
entre si todos os possíveis resultados da medição do observável (E, P, X;). 

Em MQ, pode-se distinguir dois tipos importantes de relações entre observáveis. 
Analisemos dois observáveis, 4 e ғ são compatíveis se as propriedades (E, Р, Х;) de todos 
os pares (4, rj — qi representa qualquer propriedade do observável q, r; qualquer 
propriedade de 7 — forem compatíveis ou de mútua exclusão. Sendo o último caso, o par de 
propriedades que não correspondem a qualquer situação fisicamente real ou previsível. 

No entanto, para que dois observáveis q e s sejam incompatíveis basta que, entre 

todos os pares de propriedades (q; sj) exista um que seja incompatível. 
Um exemplo de incompatibilidade é o momento linear (P) e a posição (xi). A energia de 
uma partícula livre (energia cinética) e a sua quantidade de movimento são também 
observáveis compatíveis. Neste caso, os pares de propriedades formados por uma 
propriedade de cada um destes observáveis ou são compatíveis ou excluem mutuamente; a 
exclusão múua aplica-se a situações fisicamente absurdas, por exemplo quando, uma das 
prorpiedades for nula e a outra diferente de zero porque a energia cinética e a quantidade de 
movimento estão relacionadas entre si através de uma dependência funcional, К = p^/(2m). 


1.7 Valor Esperado 


Se considerarmos, uma partícula e onda associada, a função V, e se essa função 
não se anula num intervalo entre r e r + dr, na medida de sua posição há uma probabilidade 
finita dessa partícula ser encontrada. Não podemos, atribuir a coordenada um valor bem 
definido, no entanto, é possível especificarmos uma posição média da partícula. 

Imaginemos a medida da posição da partícula no instante t, a probabilidade de 
encontrá-la entre r е r + dr é dada pela equação. 


ТЕ ЖЕЛІЛІ [16] 


onde P,.; é a probabilidade de encontramos a partícula. Repetindo essa experiência a uma 
certa frequência no mesmo instante e registrando os valores de Р), podemos usar a média 
dos valores observados para caracterizar a posição da partícula no instante t. 

Este valor é representado por <r>, valor esperado da coordenada ғ. Veja abaixo 
como podemos demonstrar matematicamente este cálculo. 


<r>= f rP, pdr [17] 


—00 


como: PF, = V cV, 
substituímos [16] em | 17] temos que: « 7 >= [rur eov, dr [18] 


—00 


1.8 Fluxo de Probabilidade 


A probabilidade de encontrar uma partícula a uma certa posição x e num tempo t é 
dado por 


Pos) = VU 


Sendo o fluxo qualquer grandeza que varia o tempo. 


ОР 0 
IT „Í = — t ed 19 
a, 09D = lY ea en] [19] 
ОР „ дү ду * 
E (xt =W eo A ТЕМЕН 20 
РЗ V ( s ы) V (xi) et (x,t) [ ] 
2 2 
0 ow 


consideramos que ela esteja livre de V(x). 


Parte real: 
. 2 
ih 0 a = ду [21] 
2m Ox Ot 
Parte imaginária: 
. 2 ж ж 
ih zi D ду [22] 
2m Ox Ot 


Substituindo [21] e [22] em [20] 


ôP ЖЕЛЕ ih дү 
оре = V (x) £ ДЕ Kaal- > | 


2m Ox? 2m ox? 


oP ol h „ дү дұ” 
,) = — x) —— Vu —— 
ot сы) ШЕН ) Ox Ven Ox | 


h ¿0 ду" 
e [ у оғ Л [23] 


Sendo 5 o fluxo de probabilidade 
ӨР os 
X,Í) = — x,t) ou 
A (x,t) EE (x,t) 


E ne (ийй [24] 
Ot Ох 


Sendo esta a equação da continuidade. 


1.9 Postulados da Mecânica Quântica da Função de Onda 


D A todo sistema corresponde uma função de onda Y (x,t). 

П) А seguinte correspondência existe entre as variáveis dinâmicas е os seguintes: 
Variáveis Operadores 
х,у,2,1 X, y, z, t 
F(x,t) F(x,t) 


P,, Ру, ou P; ho h д ои? 0 
їдхїду 10% 
Е h д 


HD V (x,t) (onde x representa x, у, z) e suas derivadas primeiras, devem ser finitas e 
unívocas. 


IV) Y (x,t) deve ser normalizada”, isto é. 


+оо +оо +оо 


BERATEN 


—00 —00 —00 


dv - dxdydz 


ou simplesmente 


* quando tem-se uma função = 1, a probabilidade é 100% de encontrar a particula numa região, isso significa 
que esta função está normalizada. 


[ечат 


—00 


ou mais simplesmente 
«WW >=1 


V) А cada observável corresponde um operador А. O valor médio <А> ou valor 
esperado deste operador e definido como. 


< А>= ШЕГУ? 
<А Mi ЕШ! 
ou 
< А>=< Pay > 


porém, se Y nào for normalizada, então <A> é definido como 


«тА > 


«A» 
< YY > 
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Capítulo 2: Equação de Schródinger. 


2.1 Introdução 


Em 1926, o físico austríaco Erwin Schródinger (1887-1961) publicou quatro 
trabalhos nos Annales de Physique Leipzig nos quais desenvolveu a sua famosa Mecânica 
Quântica Ondulatória, cujo resultado principal é a equação para as órbitas estacionárias dos 
elétrons atômicos, a igualmente famosa equação de Schródinger: 


л? 
Ее =0 0 


Em relação ao trabalho de Bohr, o trabalho de Schródinger foi bem mais completo. 
Uma vez que prevê também o seguinte: 


е As autofunções são correspondentes a cada autovalor. 

e Prevêo cálculo da probabilidade de um determinado estado. 

e Prevé o cálculo da probabilidade de transição de um estado para outro. 
e Calcula os momentos angulares orbitais. 


A equação de Schródinger nada mais é que uma equação diferencial de segunda 
ordem,a qual podemos aplicar para um sistema como o átomo de ІН e calcularmos os seus 
níveis de energias correspondentes. Historicamente foi o primeiro sistema que Schródinger 
tratou, onde os autovalores de energia são os mesmos que previstos por Bohr. 


2.2 Construção da Equação de Schródinger 


Primeiramente temos que resolver um problema de dois corpos. Neste caso, 
podemos reduzir o sistema de dois corpos a um sistema de um corpo, considerando a 
massa reduzida do sistema: 

__тт, 


[2] 


m, + m, 


Este é o termo que é introduzido na equação Schródinger e podemos adquiri-lo através do 
cálculo do centro de massa e para acharmos o centro de massa temos que fazer a seguinte 
igualdade myr;=myr, е se variarmos г logo teremos o centro de massa. Considerando 
m>>m; 


ANS [3] 


r=n+h [4] 


h=r-h5;h=r-n [4.1] 
Substituindo [4.1] em [3] obtemos, 
n- т; ren= т , [5] 
nm, + т ті +m 


Como o próton e o elétron estão translacionando е girando com velocidades próprias. 
Podemos obter a expressão da energia cinética total do sistema е a velocidade angular. 


2 2 
_ my У 


Ep = + 6 
xk ==, 3 [6] 
у= g ei 
n 5 
1 


1 
sendo I- (ти? + mar?) [8] 


Substituindo [5] em [8] 


j-l| тт |2 [9] 
2( т, +m, 
т, + mM, 
1 
Ey = = што? [11] 


r 1 r r : r . 
No átomo de Н o núcleo é massivo e o elétron tem massa reduzida y dada pela 
equação [10] de tal maneira que gira em torno do núcleo estacionário. 


elétron 


Figura 1: representação do átomo de hidrogênio 


O potencial de interação (coulombiano) do elétron-próton é dado pela relação 


Ze 
ос ses my = 12 
Але,(х? + y’ + gu 12 


(ул) 


Onde е = carga do elétron, Z = carga do núcleo ( para о 'H, Z = 1) 


А eq. de Schroedinger dependente do tempo para este sistema é a seguinte relação, 
conforme mostrado abaixo: 


wE xot 2 x OE oso 
ШЕ ! y? | 22 КЕЗЕ) Hox, y2) х,у) =iħ РУ 13] 
ou numa notação mais compacta: onde se consideramos о operador laplaciano 
О @ Ə 
2 
Жш E [14] 
X y 2 
Assim a eq. de Schroedinger pode ser reescrita assim: 
л? os ov 
— VP + VY = ih — [15] 
2u ot 
Cuja solução é dada: 
-iEt/ħ 
Payiz Фоул) [16] 


Onde yy; é a função onda independente do tempo dada pela resolução da eq. de 


Schroedinger independente do tempo a seguir: 


Ro, 
с pd Gon “Роуз Payo Еа) [17] 
ou 
nio ә? 0 
= ES Ж yl + po: pu + Vs қағаз = EN (х.у,2) [18] 


Para resolvermos a equação acima é melhor escreve-la em coordenadas esféricas e usar о 
método de separação de variáveis 


=... O elétron 
I 


Figura 2: representação da interação coulómbiana entre о 
próton e o elétron, aplicando as notações de coordenadas 
esféricas para resolver a equação de Schródinger. 

De acordo com a Figura 2, r é o raio vetor posição do elétron, 


r= ух + y? + z? [19] 


Ө é o angulo polar, o qual cresce a partir do eixo z para o plano xy 


Ө = arccos e [20] 
yx + y2 +2” 
e ф é o ângulo azimutal, o qual cresce de x para y. 
Фф = arctg > [21] 
ES 
х —rsenÓOcosó 
onde temos y = rsenGsenó 
z = гсоѕ0 
E о laplaciano em coordenadas esféricas é dado assim: 
1 2 
quu. A д no y a = [22] 
r Or Or r веп0 00 00 r sen 009 


A equação de Schroedinger em coordenadas esféricas fica assim após aplicarmos o 
operador laplaciano esférico: eq. [18] 


à |10(.,0v 1 ð ov 1 07% 
+ sen Ө + tV, MS = EY.. 
2u Е дғ Ë or | r? зеп @ A 2 r? sen? Ө дф? | шала C) 


[23] 


Nota importante: Podemos escrever a função de onda Y, „у Como um produto de três 


funções tal que 
Vlr, 9,0)= К(ғ)Ф(ф)6(0) 


1. A primeira R(r) = dependência em г ( radial afastamento do elétron) 
2. A segunda Ө(Ө) = dependência em Ө (mostra a posição polar do elétron) 
3. A outra função Ф(ф) = mede a posição azimutal do elétron na sua trajetória 


Podemos realizar esta representação porque há um campo de força central na equação. 


Aplicando agora a separação de variáveis. A equação [23] pode ser reescrita (veja notação 
da derivada e o que deriva parte constante fica fora da derivada) 


тераа, IR 2l 249) IRO do 


r + sen + +V. ROD =EROO 
2u| r dr dr) r?seng do dO) ү?веп?0 аф? | Va 


EF 2 
h 99 4 (ad), IRO dls Du IRO о | riomoo- Enos 


2u| г? dr ағ) r °sen0 do do) ү?веп?0 dp? 
24] 
2 
Substituindo o potencial V(r )=- 1 E e divindo por ROD a equação [23] fica: 
ле,” 
2 2 2 2 
sen*0 d [7 A em 0 d Ёз 424 2 en! 0 pE __ 1 а“Ф ЖЕТЕ; 
R dl dr Ө do do)| m Am&yr Ф dy? 
[25] 
2 
„кн u --т? (constante) 
Ф dy 
2 2 2 
Иа ғ a сш, |1099, p D Ey =-m? 
R drN dr © do dO h° 4zeor 
[26] 


Multiplicando por r^sen? 0 todos os membros fica 


2 2 2 2 
sen” 01 36 ба [seno o) + 1 E +28 2 sen? 9 Ey E Si 
R dr dr Ө do do) Ф do h 


Ou 


мы кд РОВ: 2452 2 2 
ғ o 0 lr zx r n 0 d CAE r m : = NU fut -0 
Rr^ dr dr) Oríseng 40 dO) Фк“веп“0 аф 


Isolando os termos em q (passando ao segundo membro) 


2 2 2 2 
беп A J 0 pena] + E isento ps __14Ф 
В drV dr e de do )| h 4xeyr Ф аф? 


Note : do lado direito а equação depende só de q е do lado esquerdo а dependência é em Ө 
e r. Como temos uma igualdade a constante de separação (-m?) deve ser a mesma para 
ambas as equações. 


2 2 2 2 
enr E Je no) eset sm 
R dr dr Ө do dO h 4z&gr Ф dy 
2 
m. =-m” (constante) 
Ф ар? 
2 2 2 
— eg б Aem ш ЕШ NM. sen? 9 E + — --m? [27] 
R dr dr Ө 40 do n? Am£gr 
2 2 
10 = ой б Sem =( cujas raizesA = tim 
Ф dy dy 
logo a solução é 
Ф(ф) = Ne“? [28] 


que é a solução particular, sendo т número quântico magnético. Ф(о) é a função onda 
átomo de ІН que contém toda a dependência em função do Potencial q. Onde q varia de 0 a 
27. Normalizando a função (o) 


(Ф|Фу-1 
2л i 
(D|D) = Гуе"? Ne“ do =1 
0 
27 


(Фф) = М? || е" ¿mo go =1 [29] 
0 


ой М'[2л]=1 => Ns 
logo a função vale 
Ф(ф) = е [30] 


como D(p) = P(g + 27) > ela é única 


ei"? =e"? e™? 1 = 927^ 1 = cos 2zm + ¡sen 27m 


veja que satisfaz se m=0+1+2+3... m = inteiro (número quântico magnético) 
escolhe a solução À = im (movimento de giro elétron sentido horário). 


Assim temos até o momento 


Y (r,0,p)= 


L RDC) [31] 
т 


Tomando a equação [27], 


2 2 š 
ш 0 z(n Ae 0 d Бо е fem 


R dr dr Ө do 40 n? 


Dividindo por sen”0 fica 
2 2 

: (е A) ша. Er] pH D pa Ca] m [32] 

Rdr dr O 40 dO h? 4zeor sen? 0 
Separando os termos da R(r ) da função O(0) fica 

2 2 

Lafe R), 2, pass ШЕ E [33] 
R dr dr h? 4zeor sen20 Ө d0 dO 


devido a igualdade a constante de separação é a mesma 


2 2 
sao exe [P E 2 ap (no j- it «n [34] 
R dr dr) h? 47E0r : 10 


2 
Id б A 221 EE dou [35] 
R dr dr) m? Az£gr 


m? als do 


— — епӨ — |= (1+1 36 
sen?g Ө 40 48) VIS [36] 


Divide а [35] por 7? e arrumando fica 


2 2 
solo São Е wanu an 
r? dr dr) m? Az&gr 


ои [37] 
2 2 

l (eE), 28 | p s +1(1+1)|Ё=0 

r? dr ағ h? Am&gr 


Esta é a equação Radial — cujas soluções são as funções de Laguerre. Apresentaremos a 
solução que é: Solução particular: 


2 3/2 Me 
n. = 2) e CE Lato) [38] 
na 
ME соу em [a DIT e" 
ases o | 2 P ЕРЕН 
быы" 
па 


n-123., 1=0,1,23., ш-1-4%1,,4-1,1 
1 3/2 
Se n=1, /=0, Rut)=(*) de 
a 


2 


Veja que se n=2, /-1 temos dois /=o el=1, Nota а = 
mq 


1 3/2 z 1 3/2 $ 
Км (r) = ЕЗ (2 ==, 7 e В, (ғ) = ЕЗ тИ 


2а АЗ 


А outra equação: 


2 
I mete) I(I -1)- Ме e=0 
sen Ó 40 dO sen? 0 


As soluções desta equação são as funções associadas de Legendre. Os harmônicos esféricos 
- são soluções das equações diferenciais: D(p).0(0) representadas por 


ү (0 ) 21 i 10 P т)! Ü P” (0) imp 
9) 3| ————— e 
BS U SED Ы 
d" p. 
onde X P/"(0)-sen L Sep. 19) 
d(cos ө)" | 


_ 1 4 een] 
MOSS "m 


onde £ = cos 0 


А solução geral para o átomo de Hidrogênio ё: У, (r,0,@)= R,(r)Y,, (0,9) өлме 


Que pode ser escrita assim: 


EVI - ы м nz D ON +1 [a+ DIT £^ -(i Hh) Egt 
Fam (r,0,9) = 2) е 2 $ Lo. o 2n[n "T b P 1)? (n 1 1 р)!01 ETÀ ppt (Ө, pje 


2r 
Com ¿== 
na 
n =1,2,3... número quântico principal 
| =0,1,2,3... número quântico de momento angular 
т =/¿-1+1,...1-1,1 número quântico magnético 


2.3 Números Quânticos 


Cada conjunto de л, /, m define uma função de onda — que é um estado eletrônico do 


átomo. 

n=1 camada k 
n=2 camada 1 
n=3 camada m 
n=4 camada n 


n=5 camada o 
As funções de onda de cada camada são chamadas de orbitais. 


a) para cada valor de n, há n-1 valor de / 


n= = 0 
п = 2 [= 
n=3 - 2 


п define o estado de energia 
b) para cada valor de 1 há (21+1) valor de m. -/< m < l 


[e т-0 


т--2 

m=-1 
[=2  -I<m<l m=0 

m=1 

m=2 

m=-=1 
[=]  -I<m<l m=0 

m=1 
[2-0 m=0 


Concluindo, quando resolvemos a equação de Schródinger para o átomo de 
hidrogênio, obtemos as energias para os níveis energéticos (En) е os orbitais atómicos 
(Paim). Temos então um diagrama de níveis energéticos no interior do átomo, i. e., 
dizemos que um átomo se constitui de um sistema energético quantizado. Na Figura 3 
demonstramos soluções da equação de Schroedinger radial. 

E foi no Annales 80, que Schródinger publicou o seu estudo do efeito Stark através de sua 
Mecânica Quântica Ondulatória (1). Um dos seus primeiros trabalhos utilizando a nova 
teoria. 


Hydrogenic radial wavetunctions - 2s orbital Hydrogenic radial wavefunctions - 2p orbital Hydrogenic radial wevclunctions - Эр orbital 


---- E E - +— | 


210 


2х10 


в - — „сэ!“ 


R 1 =) demís 2/4) bee се22і ag | 1 (2\5 1 E 
Rural] 9-2 2егіс 2/4) where p= 22012, А 1 E KA К," — || -Aep 2/6) where o= 22r i ag 


i perp- old) where ә” 22r las 26) 


06) (а, 


221 


14 


Hydrogenic radial wavefunctions - 3d orbital Hydrogenic radial wavefunctions - 1s orbital Hydrogenic radial wavefunctions - 3s orbital 


3 


212. АРТЫ 


x 
Ra” aix] (д exp(- ol 6) where om 227! а, X 
8130) (а, 2 
| 8.42) ехр(- ol 2) where ,0=2Zr1 ag NES 
2-1 а, зх 
2-1 
2-1 


Figura 3: Soluções da equação de Schroedinger radial. 


2.4 Função de Onda na Concepção de Max Born 


Bom (1928) deu um passo a esta dificuldade propondo uma interpretação estatística 
das funções de ondas do elétron, à qual, devido às inúmeras vantagens apresentadas, tem 
sido amplamente aceita. Born supôs que as ondas não têm existência real, e assim, as define 
como ondas de probabilidade. 


O produto FP ou | em um ponto representa a densidade de probabilidade de 


encontrar o elétron, ou um outro corpúsculo qualquer, em um ponto x, y, z, num dado 
instante / e igualdade 


¡Y dxdydz = ЕЛ 


Representando a densidade de probabilidade de encontrar о mesmo elétron em um 
elemento de volume dv, e também o número de elétrons dentro do mesmo volume. Esta 
interpretação teve um pleno acordo com as condições de Schródinger (5). 

Em processos vibratórios o conhecimento da amplitude é importante como o 
conhecimento da frequência própria; analogamente, é de se esperar que, em mecânica 
ondulatória, esteja ligado um importante significado físico à função de onda V ou antes, ao 
quadrado do seu módulo, visto ser evidente que o valor instantâneo da própria função 
oscilatória nào pode desempenhar qualquer papel em virtude da sua alta freqüéncia. O 
motivo por que se toma o quadrado do módulo é que a própria função de onda (devido ao 
coeficiente imaginário da derivada em ordem ao tempo da equação diferencial) é uma 
quantidade complexa, enquanto as grandezas suscetíveis de interpretação física devem 
evidentemente ser reais (6). 

Suponhamos que no estado caracterizado pela função de onda V, se efetua uma 
medição que conduz com certeza a um determinado resultado, e que o mesmo fazendo o 


2 % 1 
5 (6-2,0+= 4) exp(- ol 6) where 2-2 
a 9 


estado no estado V», conduz ao resultado 2. Admite-se então a combinação linear de Y; e 
Y), о que significa que toda função de forma CY; + CY) (C, e C3 , constantes) representa 
um estado em que a mesma medição pode dar um resultado 1 ou o resultado 2. Podendo 
afirmar, que se conhecemos a dependéncia dos estados com respeito ao tempo, dependéncia 
a qual é dada pela função VV; (x,t) e em outro, por Y2 (x,t), pode-se notar que qualquer 
combinação linear destas dá também a possível dependência de um estado do tempo. Estas 
afirmações constituem o conteúdo do princípio de superposição dos estados — um princípio 
positivo fundamental de mecânica quântica (7). 


2.5 Valor Esperado da Função de Onda 


Consideramos uma partícula е onda associada, a função “(с e se essa função não se 
anula num intervalo entre r e r + dr, na medida de sua posição há uma probabilidade finita 
dessa partícula ser encontrada. Não podemos, atribuir a coordenada um valor bem definido, 
no entanto, é possível especificarmos uma posição média da partícula. 

Imaginemos a medida da posição da partícula no instante t, a probabilidade de 
encontrá-la entre r е r + dr é dada pela equação. 


F^ V. (VC) [16] 


onde Р.) é a probabilidade de encontramos a partícula. Repetindo essa experiência a uma 
certa frequência no mesmo instante e registrando os valores de P4, podemos usar a média 
dos valores observados para caracterizar a posição da partícula no instante t. 

Este valor é representado por <r>, valor esperado da coordenada ғ. Veja abaixo 
como podemos demonstrar matematicamente este cálculo. 


+оо 
<r >= f РР, dr [17] 
como: P, = V (aW, po 
+00 
substituímos (16) em (17) temos que: < 7 >= [rr cwn dr [18] 


—00 


2.6 Aplicação da Equação de Schródinger 


Tabela 1: Exemplos do operador Hamiltoniano para o movimento de uma 
partícula de massa m em diferentes campos de força definidos pela função 
(operador) potencial V. 


Operador 
a) Partícula livre e— E 2 42 
m dx 
b) barreira de A 
( ) | V=V, Н- d (х<0; х>0) 
potencial Б 2та 
» 2 2 
0 а u Hi== pe, Vo (O<x<a) 
2m dx 
1 2 2 
o HM Oscilador K у-(1/2)Кх? Н- ла 1 Kx? 
armonico 2m dx? 2 


Tabela 2: Exemplos do operador para átomos e moléculas. 
Operador H 
(a) Átomos de |^... Ze 
um elétron шеттей Er 


(b) Átomos de л? s We e 
muitos elétrons mp 2 uS 
E FE 


(c) moléculas RO ees 2, ZaZa 
ка че. XJ a DD - po E 
i= i j»i ti а В>а 


2.6.1 Poço de Potencial não relativístico 
Considerando um potencial degrau unidimensional definido por 


V(x)=0,x<0 
V(x)=V,,x > 0 


Supor que uma energia incidente da esquerda para a direita tem energia E = 4Vo. Neste 
exemplo podemos calcular a probabilidade de que a onda será refletida (coeficiente de 
reflexão). Analisando o contorno do problema, para V(x) = 0 quando х<0, o potencial na 


região I é nulo e consideramos a partícula livre. A equação de Schródinger se reduz à. 


o, 


ox? 


+k? EY, =0 


As raízes da equação: 
Д2 +k? =0, X = Е, A=tik, 


Y, = Ае Ms Ве cH 


Ae" — onda incidente 
sendo 


Be 8” sonda refletida 


Solução geral 


v, =P e Eh 
Меке (4e + Be tt" егіні 
W,=4 еїйїх-Ё!л) +В (txt Er/h) 


Na região II há um potencial V(x), sendo assim a partícula sofre uma ação do potencial. 


2 
m 
ho 
cmi Yy = [E - Vo y 
m 
2m 
V^, Pa Voy =0 
2m 
Considerando kj = zz E А 
VP, 42%) =0 
Achando as raízes da equação 


Д2 +62 =0, 4 =-k2, A=tik, 


vw, =(Ce2 + De) 
A solução aceitável é 
Y, Es Coe 


Solução geral 
m iEt/ h 
Py = Фе 
Es Жох | iEt/h 
Ур = Се “е 


Y, = Се бох-Е/а) 
logo temos as duas funções para as duas regiões I e II 


Y, = Ае "1% + Bei” para x «0 


Y, = Ce"? para x 20 
Analisando as condições de contorno 


ik —ik jk: 
a) YP; |0 Yy lico Ае!“ + Be  =Ce > 


А+В=С deve ser contínua no ponto х=0 
b) A sua derivada também 


a T ax 


ik Ae" — ik Be “1 = ik, Ce"? 


ik (4- B)= ik,C 


(4-B)- Ë c 
kı 
Montando o sistema 
Achando B Achando C 
A+B=C A+B=C сазы, вас 
"шый ma Matan Es n 
k, k, аа 
24=-2с+С 2-1-8) C 2\ k 
! l k. +k 
С=-—2В|——- 3 
2 kı 2 К, 


Substituindo os valores de A, B e C em Y, e Yy 


y, = Spem LEE зі х<0 
1 


Y, = C|| k +k, ¿xr Im) y k-k eia Eun) | x<0 
2 k. k. 


фр =-2B k-k еї(йзх-Е!!/л) 
L К 


Obedecendo ao princípio da complementaridade temos a combinação linear do dois estados 
da função de onda 


Prota = Y; + Wy 


Хт = с ktk, ейбах-Е/)) ЖА k-ky eg Vax-Et/n) -2B ki _ k, gas Ett) 
2 kı ki 


kh 
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Capítulo 3: Método de Hiickel. 


3.1 Método Semi-Empírico 


A equação de Schródinger para moléculas não pode ser exatamente resolvida. Isto 
ocorre porque a equação exata não é separável. Uma estratégia é de assumirmos que 
podemos permitir de escrever uma forma aproximativa de equação de Schródinger para 
moléculas e principalmente que esta nova equação de Schrôdinger se torne separável. E há 
então de procurar como resolveremos esta equação separável. O método de Hiickel é uma 
possibilidade e de campo autoconsistente é um outro. 


3.2 Elaboração 


Separação o — 7 


Para o tratamento de moléculas com elétrons т é possível escrever o Hamiltoniano total 
como: 


HeH,-H,-H, [1] 
o método de Hückel extendido é 

Н(о)%(о))- El) Y (o) [2] 
e o método de Hückel simples é 

Hr) (2) = E| Y) [3] 


O Hamiltoniano que descreve a energia do elétron л é 
1:22. 
иче ү (D +/, (0) [4] 
Sendo о Hamiltoniano total 


HI"q)-Y H,O [5] 


HRS 04,0 [6] 


A função Y, é antissimétrica de n órbitas moleculares л. А função antissimétrica 
Y, descreve partículas que obedecem a estatística de Fermi-Dirac, para Férmions. O 
método de Hiickel é um método classificado como semiempírico. Este pode ser 
simplificado no máximo a resolução das equações (E) = [v mac , Ou seja ele utiliza do 


método variacional para solucionar problemas de moléculas com orbitais л, e para isso ele 
faz algumas considerações. 


3.3 Considerações do Método de Húckel 


i) S; = д, onde 6 é o delta de Krónecker. 
ii) H; =œ onde o, é um parâmetro coulombiano 


iii) A, = B sei е pertencem a átomos ligados entre si е 
Estas considerações são realizadas sobre a determinante secular 


|н, - ES; 


=0 [7] 


3.4 Aplicação do Método de Húckel no Etileno 


H52C——CH5 
А determinante secular Н, para o etileno Y, = 2р, e Y, =2p.(2, será 2x2 


Hy - ES, Hy - ES, 
Ho - ES»; Н, - ES» 


Utilizando-se as condições de Hückel 


а-Е p 0 
P @-Е| i 
В 
а – Е 1 
B go 
з HE 8 
ГА 
x 1 
-0 
1 x 


Esta é conhecida como determinante de Hiickel. 
x?-120, х= +1 


рагах=1 E = а – В ерагах=-1 Е = а + f 


» 01702 


Quando normalizamos as funções [ча т =1 nós achamos os valores de сі 
ж ж 
flet, +c Y, Jo, +e Y, kr =1 


ci | Widr + c; | Wi ade + er Y) dr + c; |247 =1 
0 0 


с = 


al- 


logo achamos a funções V e Ҹә 


Y, 20.707 (p) +0.707%(р.,) 


V, 20.7079 (p) — 0.707 (p.,) 


34 ө . ө 
с Etileno 
24 
1- 
0- 
ө ө 
ө ө 
14 ...... 
Т Т Т Т T T 
2 1 0 1 2 


Figura 1: comportamento da função de segundo grau do Etileno 


3.5 Aplicacáo do Método de Hückel no Butadieno 


A molécula 1,3 Butadieno é plana e tem preferência como isômero trans е tem 
fórmula 


Нүр -ES Ho -ESj; Hj3-ESj3 Hi4-ESi4 
Нз — ES» Hoy — ES» Нҙҙ- Еәз Нэд – ES @ 
Нзү-—Е$зү Нҙ)-Е5ҙ; Нзз – Е5зз H34 — ES34 
Ha —ESg На -Е5)) Наз – Ез Над — ES44 


Fazendo as considerações do método de Hückel 


а-Е p 0 0 
p «-E 0 E 
0 p «-E 
0 0 p «-E 


1 0 0 
В 
p EE. q 0 
á а-Е o 
0 1 1 
p 
0 0 еее? 
В 
TE 
x100 
lx 1 0 
-0 
0 1 x ] 
001 x 


quando se tem uma determinante diagonal simétrica, utilizamos a fórmula abaixo 


IT a дс)? 023) 


j=! 


sendo j = 1,2,..n 
x, l4, 0 
l. x 1 0 
с а b =0 
0 l x 51 
0 0 L. x, 
4x* - 3x? +1=0 
as raízes são x = 40.62 e х= +1.62 
E, =а + 0.620 
E, =а — 0.628 


E, =a –1.628 


4n 


a — 1.62 В 

31 
a — 0.62 В 

С,, 

2n a + 0.62 В 

Іл 
а + 0.62 В 
x 10 0\ (а \ (0 
1 X 1 0 C5 Шы 
01 x I[le| |0 
00 1 х) (с) lo 


су + Xc, = Ü 


Agora que temos o sistema montado para o Butadieno, podemos calcular todas as funçóes 
de onda relacionadas com as ligações do Butadieno. Iniciando a calcular V^, o orbital de 
menor energia, utilizamos x= -1.62, vem: 


C5 =C3 


с=с 
Agora com a condição de normalização de “Р, (para S;=0) 
[CET +07 Pa +с;Ф; + e V, 0, V, +e, Y, +c, Y, dr =l 


c2 63 


c? +e? +e? +c =1 
c? +(1.62c) + (1.610) + c2 =1 
1240" =1 5050.72 
substituindo c em V^, temos o resultado da primeira combinação 
Y, =0.372y, +0.602y, + 0.602; + 0.372y, 


e assim para acharmos Y,, Y, e V, devemos substituir no sistema acima apresentado. 


“0 Representação 
de 4x*-3x+1 


Coeficiente C 


Figura 3: comportamento do determinante do butadieno 


3.6 Aplicaçáo do Método de Hückel para o Ciclobutadieno 


A molécula do Ciclobutadieno difere da molécula do Butadieno porque esta possui 
os carbonos ligados 2 em 2. lembrando que a determinante secular é 


Hj -ESj Hj? -ESjj Ніз – Ебі 
Ho —ES Hn —-ES» Hz — ES 
H31-ES31 Нҙ)-Еӛҙ) H33- E83; 
Ha =ES4 Ha – Ед Наз -ES4 


а-Е 
B 
а-Е p 0 В 
p а-Е 0 Ж | 
0 B а-Е B E 0 
P 0 В а-Е1 
1 
x 1 0 | 
а-Е lx 1 0 
sendo х- - 
0 1 x 1 
101 x 
x=0 е х-ЖЮ 
Ез =а a 
onde E; е Ez são degenerados 
E, =а 


о-2р ==> 


о--2р 9% 


x 1 0 l| íc 0 XC +C, + C, = 0 
1 x 1 Offce] JO сү + хс) + c4 = 0 
0 1 x 1| с; 2/0 €, + хез + c4, = 0 
1 0 1 x) Ке 0 с “су + Xc4, = 0 


Figura 5: comportamento do determinante do ciclobutadieno 
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Capítulo 4: Método Variacional. 


4.1 Introdução 


O método variacional trata-se de um princípio importante em Química Quântica, 
não só porque permite obter soluções aproximadas para a equação de Schródinger como 
também porque fornece um critério, o valor mínimo de energia. Considerando a equação de 
Schródinger. 


НҮ = EY [0] 


onde Y e E são autofunção e autovalor do Hamiltoniano, são soluções exatas desta 
equação. Normalmente o Hamiltoniano é de tal forma que a equação [0] não é separável, 
qualquer que seja o sistema de coordenadas. Por essa razão não é possível pensarmos em 
obter uma expressão analítica exata para V. Passemos então a examinar as propriedades 
da funcional E(w), definida como 


i | Ф“Н%ат 
f par 


(E 


s[wmrac [1] 


Esta é a integral variacional. Uma função qualquer, que seja normalizada, bem comportada 
e que satisfaz as condições de contorno do problema, podemos mostrar que esta função 
admite que o valor esperado do operador Hamiltoniano para uma função de onda 
aproximada é sempre maior que a energia exata E, correspondente a esse Hamiltoniano. 


f Ф“НЧат>Е, [2] 


Em outras palavras, a função de onda aproximada conduz à menor energia. Assim sendo, os 
coeficientes c, nas orbitais moleculares V^ são determinados de modo a conduzirem ao 


valor mínimo de energia E associada com essa orbital molecular. Em relação à obtenção de 
soluções aproximadas, o que se faz na prática consiste em partir de uma função aceitável, 
com vários parâmetros ajustáveis, e variar esses parâmetros de modo minimizar E pois, 
segundo o princípio variacional, conseguirá obter-se uma melhor aproximação para a 
energia do estado fundamental Æo 


[De HD par yxa Í | we AY, Yee H yy 
_ k k _ k j _ Ki 
ГУ се У ay УУ скс | PF; Y eie;Sy 
k k k j kj 


(E 


Provando a afirmação [2], 


=> ra s, [3] 
k 


Substituindo [3] no lado esquerdo de [2] 


| v'mvac- [Y a W HD a V dc- [Уа V, у a BP dc 
k j k j 


(Ғнеа-Гуауж у aE Y dr= у у ауа,Е,| V, V ,dr - У у a;a; E, V ó, 
k j k j k j 


[Y mrac- Y у aja E,V ó, [4] 
k j 
, | 1 ве К-/ 
onde б, é o delta de Kroenicker, 
0 se Кж 


. ж 
como temos ó,; podemos conseguir o termo EN aja, E, 
k 


ж ж 2 
k k 
sendo que o nível k é sempre maior E, > Ер, onde Eo é o estado fundamental do sistema. 


[par - Уаа,Е, ЭЛЕ -E Yl [6] 
k k k 


Através da normalização 
[ V'HYdr -1 [7] 


= Жаку ат Y vaa; | Wy de 
k j k j 


agora substituindo W 


УУ аа, Y usai [8] 
£ j 

abrindo a somatória em / 
Y Y aia =D lal [9] 
k j k 


Substituindo [9] em [6] 


[v Hear =E Y Y ака бу [10] 
k j 


pe ass Е, [11] 


Isso mostra a prova. Por isso temos que ajustar bem a solução de 'Y muito bem para 
atingirmos o estado E, (valor real) que é o estado fundamental. Ajustar a função Y é 


ajustar a função de base computacionalmente. 


4.2 Extensão do método Variacional 


O método variacional apresentado na seção passada tem duas limitações. Primeira, 
ele provê informação de energia somente sobre o estado fundamental e da função de onda. 
Segundo, ele provê um estado acima do estado fundamental somente. Nesta seção 
gostaríamos de estender o método variacional. 

Considerando que nós desejamos estender o método variacional para obter um valor 
de energia estimada para o primeiro estado excitado. 

Demonstrado que a função variacional normalizada [5] e [8] 


| puq. S ja ne e jap ara ХЫ =1 


к=! ka 
onde a, é a expansão do coeficiente em Y = У a, V, . Tendo a, = (Ф, |У). A função Y 


k 
que é ortogonal para o estado fundamental real da função de onda Фо. Então nós temos 


a, = (9, |W)-0 е 
e 
(ғнғаг- Y s E e [vide Уа, =1 [12] 


k-2 k-2 


рага k>1, temos E, 2 Ее la, E, 2 la, | E, 


Y la, E > Ma E = E Ya = ЕЁ, [13] 
k=1 Еа k=1 
Combinando [12] e [13], temos o resultado: 
Гғ“Нғаг>Е, ‚ ве [#°ЙФат=0 е Гғ“Нғаг-1 [14] 


a extensão do método variacional para altos estados excitados acima de Ej. 


[v Arva: i | | 
— C у 2 Ep SO 0 T= 1 t=""= k т = 
[rvar >E [rodar - | vid [rimar =0 [15] 


4.3 Função Variacional Linear 


Um tipo especial de função variacional é usado largamente no estudo de moléculas, 
é a Função Variacional Linear ou Método Variacional de Raleigh-Ritz. Uma função 
variacional linear é uma combinação de n funções linearmente independente fi, f,....f, : 


Pedi tel fare Del, [16] 
/-і 


onde Y é a função trial variacional e os coeficientes c; são parámetros para ser 
determinados pela minimização da integral variacional. A função f j deve satisfazer a 
condição de contorno do problema. Nós apresentamos que V^, c; e f; são reais. 


Aplicando o teorema variacional para a função variacional linear, nós temos 
[v Wars wer Усола = у Усс, |7; Аат [17] 
j=l k=1 j=l k=l 
Agora devemos definir a integral de superposição S como 


Sa=[S; fidt [18] 


sendo jk os elementos de matriz. Substituindo [18] em [17] temos 


n n 
[wwac- у Das д [20] 
j=l k=l 
Note que 5 „ não é necessariamente igual a 5 „, desde não há razão para supor que a 


função f, são mutuamente ortogonais. O numerador da integral variacional é 


[авта y ef HY d fades» усе | unas 
j=l k-l 


j=l k=l 
e usando a abreviação 


H jy = | f; Нат [21] 


Reescrevendo о numerador 
[WHydr=5 Уссат [22] 
j=l k=l 
Substituindo [22] e [20] na integral variacional W, 


[par LIC y 


ШЕРТ 


[тєр ы oo 
[ous DDS а 


j=l К=1 


[23] 


Arranjando [23] 


Wo» уе бы Хауа а [24] 


j=l k=l per 


agora devemos minimizar W рага aproximar tão perto de Eo (W > Eg). А integral 


variacional W é uma função de п variáveis independentes c,,c,,...c, 
W=W(c,,C,,.-6,) 


Uma condição necessária para minimizar uma função W de muitas variáveis é que a 
sua derivada parcial a respeito a cada uma de sua variável seja zero no ponto mínimo. 


=0, i=1,2,..n [25] 


agora nos diferenciamos [24] parcialmente em cada c; para obter n equações: 


Y, ООК y Se; Ch 5а УУ санд [26] 


дс; 5 kl OC; 5 a Ci jsi kai 
i= 1,2,....n 
Agora 


n n ec. 
— с.с с.с а= os mno! S. 
25337 5 je ге Is S u k jk 
1 j 1 j j 


0€; 53 a j=l k= 


os cj são variáveis independentes, e portanto 


—.20 se / zi — =] se j=i 


Nós temos então 


n n 
YY, jC. jk =P Ye óS jk 5221 Ó$ д = Sm DAT 
k=1 j=l 


дс, іа j=l k=1 j=l k=1 


0 n 
2020 CrS = AS ФУ 27] 
j 


с; 


Onde nós temos desenvolvido 
S. =S =. [28] 


Resolvendo a equação [27] a partir desta relação, conseguimos compacta-la: 


556 195 з аза “Хе jy Dersa PX ik 


дс, "ci 


— n n 
eje E ik [29] 


дс, іа 


onde o fato de j ser uma variável aleatória utilizada е recolocando 5, por H , em cada 


uma destas manipulações, nós temos o numerador da equação vai 


0 n n n 
LLH к 22 еН, [30] 
Ci j=l k=l k=1 
o resultado depende de 
Н,-Ну;-Н; 


Porque Н é Hermitiano, agora substituindo [25], [29] е [30] em [26] 


W So Урнаа 
k=l k=l 
Rearranjando 


ou 
la Sue ]=0, і-і2,.ə [11] 


k=1 


Sendo esta conhecida como Sistema de Equações Seculares. A equação [31] é um 
conjunto de n equações simultáneas, lineares е homogéneas em и de ci, C2,...Cn 
desconhecidos. E pode ser representada como D;, (W ) = ІН ik — WS; | =0, onde m é a linha 


da determinante. 


Uma solução trivial das equações seculares é todos os coeficientes с, serem nulos. 


Qualquer solução fisicamente relevante exige que seja nulo o chamado determinante 
secular, constituído pelos elementos Н. — ES; 


[Hi – ES; |= 0 [32] 
ou 

D(E)=|H; ES ¡y | =0 
da equação [31] podemos obter os valores reais de W =W 1 .W.n2 ,.... нт OS quais 
correspondem Y = У, ¥ m2,- Р: Onde о valor esperado da energia Y é facilmente 
visto por 


[v mvac C^ HC 
ETE CREE 
[v Фат C'SC 


E 


4.4 Aplicação do Método Variacional 


Sabemos que agora temos a parte secular da integral variacional e podemos a partir 
daqui aplicar o método variacional no estudo de energias de orbitais atómicos e 
moleculares. 


E = | Y" par > Eo Equação Secular > Hg — ES | -0 


Para átomo as funções de base são os elétrons e para molécula as funções de base 
são os orbitais moleculares. Ex: para o Benzeno há 36 orbitais, ou seja 36 combinações 
lineares. Logo Y, será 


36 
Ж; = Уе, 
¡=1 


e a sua determinante será de 36x36. há 42 elétrons no benzeno e 21 orbitais estáo ocupados 
e outros 21 são virtuais. Sendo que o último orbital ocupado é chamado de HOMO 
(Highest Occupied molecular Orbital) e o primeiro orbital não ocupado é o LUMO (Lowest 
Unccupied Molecular Orbital). Se soubermos a diferença de energia LUMO e HOMO, 
podemos achar a energia de ionização da molécula. 

A energia total dos orbitais ocupados é da seguinte forma 


n 
Ет =D (E, +Ec+ En)? 
i=l 
Onde n é o número de orbitais ocupados, 2 o par de elétrons e E; número de orbitais, Ec 
repulsão coulombiana e Er, é a energia de troca. 


4.5 Método Variacional para o Átomo de Hidrogênio 


A partir da função tentativa Y, =e ” para o estado fundamental do átomo de 
Hidrogênio, calcule o valor de œ que minimiza a energia. 
Pa a Y 
Е= ( а | | a) [1] 
(Pa |Ya) 
Sendo o Hamiltoniano H, 


2 
ga Г (ro) 1 [seno |+ | Ó |7 
21,2 Or Or] r?seng Or Or) r?sen?Q 09? P 


ERSA -оғ n e n! 
dr = ғ? sen біғаб%ф f e "y" dy = A (a>0) 
Como o átomo de Hidrogênio depende somente de r, o operador se reduz à 
"I 1.0 ,2 0 Z [3] 
2| p2 ôr Оғ r 
viz Шз và. | 4,4 ] 22. 
r2 Or or r? 2 d\dr 
2 
Re pa LER. [4] 
rdr dy 
Substituindo [4] em [3] 
а уйш” [5] 
2 r 
Reduzindo a equação [1] 
E =(Po |H|Ya) [6] 


e aplicando o operador em [6] 


1 1 
E=-7(Ya |v^|v.)- zr; МЫЗ 


1 П E 1 
Е=- (Р |7 d E |а) - 20% | a) 
га 1 d? 1 
зш ат. | Жу Е а 


para que a variação funcional, E (Yg ), seja mínima > óE(Y,) = ó(H) =0 


d d | -ar —ar а? 2 
sendo que |9) = e )=-аје )--а Ya) e que pol x ЕР 
r 


dr dr 
teremos: 
1 
Eco 19) (Fa Pa) Z(v, |-| Po) 
2 
pote” Z 
r 2 r 
2 
Е-(а- 21) = — [8] 
r 2 
Normalizando V, 
(p |“ ed 
lembrando que [errar = т. (а > 0) 
Ге 2ar ,,2 dr =1 
222% 13 N - даз [9] 
2(4a?) 


calculando o valor médio de u 
" 


(Ша [10] 
т 


Substituindo [10] em [8] 


2 
[04 
E = (а — Z)a - — 
( ) 2 
2 2 
=й? ыи == == [11] 
ZO 


como ÓE(Y,)=0 (Н ) = 0, o mínimo de energia será obtido se derivarmos E em relação a 


о e igualando-se a zero, 


dE 
da 


-2--2-0-а-2 [12] 


Substituindo [12] em [11], e lembrando que Z=1 para o átomo de hidrogênio. 


2 
E- О qua 


Е---аиа. [13] 


4.6 Método Variacional para o Átomo de Hélio 


Aplicando o método variacional para o átomo de Hélio e encontrando a energia do 
estado fundamental a partir da função tentativa. 


3 3 
Is() = | e e 1s(2) = | —e 72 [1] 
A T 


Um Hamiltoniano de N elétrons e M núcleos é 


А 14 M | ММ? N N | M MZ 7 

НЕЕ SV =b VE LL УХ o D 
251 ka2M, EAE id ^ kk Кы 
— A / A 


220- Dp repulsáo-nn 
repulsáo-e n  repulsáo-e e Р 


No átomo de Hélio, consideremos os núcleos fixos, podemos omitir o termo de 
energia cinética dos núcleos para obter o Hamiltoniano puramente eletrônico, separando-o 
do nuclear. Sendo que há 5 termos no Hamiltoniano do Hélio 


ü --ly-vj п oda E [з] 


eletrônico 
Фар = O espacial x Ф spin 


Ҹар -1s()1s(2)xa(1)a(2), há 4 possibilidades de função de spin desde que 


negligenciamos o efeito magnético. 


A) Pap =Iso(Dlso(2) 
B) Pap =15/3(1)153(2) 
С) Pap = 189(1)158(2) 
D) Yap = 159(1)150(2) 


Onde А) еВ) são simétricas е С) ер) são antissimétricas. 


C) Yap = 1sa(D1s8(2) —IsPMIsa(2) > Yap = 


J2! 


D) Pap -1s/(DlsaQ2) – 1se()1s8(2) > Yap = e 


A 


Ise(D 1sf(1) й 
189(2) 15802) 


Isp) 180(1) si 
sE) 1sa(2) 


Utilizando C) 


Yap = 159(1)1568(2) - 152(10150(2) 
Yap = (Ls(DIs(2) L2) - (1501)15(2))8(1)0(2) 
Ya = 18(1)15(2)[0(1) 802) – 801)0(2)] [4] 


A parte espacial é simétrica e lembrando que а função de spin ја está normalizada 


Ep = BOBO- [e(1) 82) - AMA] [5] 
espacial ` - 
spin 


Verificando a ortonormalização da função de spin 


2» a AA AS 
7X leo У, POE “0505 PR I6] 
1 2 1 2 


X (ЖҮЛ УЗ PODE BOPE kP =1 
1 2 1 2 


1/2 2 1/2 2 
Sem =1 e > Mm) =1 [7] 
тұ--1/2 тұ--1/2 


Desde que as funções ое p correspondem a diferentes autovalores do operador Hermitiano 
5, , eles são ortogonais. 
1/2 


Уа (m,)B(m,) «0 [8] 


тұ--1/2 
Aplicando estas notações па verificação de ortonormalização 


1 1 1 1 
¿0:0-,0-0+,0-0)+,0-0)=1 
1=1 [10] 


A integral variacional é 


ж 


[PEE fS GODE eD- Ba) HD) -Lad 
TE 2 2 


Neste caso H não atua na função de spin e a integral variacional torna-se 


[PRY = [| [GOLDS HGD avia Y Y ad) BO) – BM [11] 
1-2 


1 
Como а função de spin é normalizada, a integral variacional se reduz a 
[wmv zi (OLD H(s(DIs(2) dv dv» [12] 


Substituindo o Hamiltoniano [3] em [12] 


[wmv = [| ШЕ [v2 -v5|-£- 2 = О [13] 


1 
"у 


ІШ + [150) 18(1)4% | so- ШІ 


Е= | 15(1) 15(2)4%› [15(1) |- 5У 


+] 1s(1) ME f 15(2) 60а +] f (15@15(2))` 2 (вуз) 
lj 


Aplicando os operadores nas fungóes espaciais 


3 3 
2 |= E de ШЕЛ 
dr; dr; | V 7 z 
3 3 
Ly) 4 je отоу | [e e; = ајә) 
dr; dr; Лл л 
e que 
qe 2 d? 2 
= | r bd —|)-q^WY 
Seem. fees 


Logo a expressão de Energia na notação de Dirac torna-se 


2 2 
1 d ру lad 1 d 1 а 


1 
а o. a | 


-z+ +z +z 


2 r; FK; FK v) 
drj i Ј ij 


E--(*| 


1 a? 1 a? 1 1 1 
=P) A] | emp) zi --zee| -mezen- p, 
rM MT 1 Ç 2 1 — A 
а а а а (5/8)a 


E zy -2га+Еа 


como óE(Y,)-ó (H ) — 0, o mínimo de energia será obtido se derivarmos E em relação a 


о e igualando-se a zero, 


UB. ador sg 
da 8 


-Xe*2--- 


3 
а--2%-- 
16 


б —1.6875 
Substituindo o valor de o na expressão de energia 


E = (1.6875)? — 2(2)(1.6875) + ža 6875) 


Е--2.8477 Hartree 


4.7 Método Variacional para Sistema de dois Spins 


PP 
ЕКІ) " 
(PI) 
Sendo o Hamiltoniano H, 
Н--уі — vol, +A, 25 [2] 


Agora analisemos a interação spin-spin e podemos ver que esta atua como uma 
perturbação no sistema de spins. E por isso consideramos esta como uma perturbação, logo 
podemos escrever o Hamiltoniano [2] da forma: 


Н-НО.н0 [3] 


onde HO =H0 + н(0) е HO = vie HO = val», o qual a sua função 


aproximada é 


Фар - Ф espacial % Ф spin [4] 


Ҹар =Is(Dls(2)xo(Do(2), há 4 possibilidades de função de spin desde que 


negligenciamos o efeito magnético. 


a) Фар =1so(Dlso(2) 
b) Yap = 15/3(1)1583(2) 
c) Фар =15&()15,3(2) 
d) Yap = 180(1)150(2) 


Onde a)e b) sáo simétricas e c) e d) sào antissimétricas. 


c) Yap = 150(1)158(2) - IsA(D1sa(2) > Yap = d 


lsa(1) 18080) =) 
21 


lsa(2) 1sB(D 


d) Yap =15P(Disa(2) - Ise(D1s8(2) > Yap = d 


ІЙ) Isa] o 
J2! 


lsB(2) lsa(2) - 


Utilizando c) 


Ҹар -1sa(D1s5Q) - 158(1)150(2) 
Ҹар = (1501)15(2))о (1) 202) — (15(1)15(2)) 400212) 
Ҹар = BODO - 201)0(2)] [5] 


A parte espacial ё simétrica e lembrando que а função de spin ја está normalizada 


1 
Yap = Om [о(1) 802) - &0aQ)] [6] 
spin 


Verificando а ortonormalização da função de spin 
1 x 1 
22 7518020) Aa] 75180002) - B(Da(2)]= 
E am? Y. Он -X DEUS Parr m 
1 2 1 2 
E B DaF PDDE BOPE kP =1 
1 2 1 2 


1/2 Р 1/2 ; 
lam, =1 e Y | Am)? 21 [8] 


тұ--1/2 тұ--1/2 


Desde que as funções a e B correspondem a diferentes autovalores do operador Hermitiano 


S, eles são ortogonais. 


1/2, 
Y a (m,)f(m,)-0 [9] 


тұ--1/2 


Aplicando estas notações па verificação de ortonormalização 
2o0-0-20-0+20-0+20-0=1 
2 2 2 2 


1-1 [10] 


A integral variacional é 


[FEES | ооо) ep) - EM] HOO) - EMA) ға 
13 J2 J2 
Neste caso H não atua na função de spin e a integral variacional torna-se 


[s mw = [ [GOLDS Hs MIC) dvd, Y Y Jew BO- BM B1] 
1292 


1 


Como а função de spin é normalizada, a integral variacional se reduz a 
[s mv = [| f (OLD н(ів(1718(2)/414; [12] 
Substituindo o Hamiltoniano [3] em [12] 
[s inv = [| BOLDY [0 +HO + HO ІШИКЕУТІР [13] 


E - | 1501504 [з Ent soam, + [10 15v 159" voi, sav; 
+ f f (15(015(2))` | Ji, À, fisco di 
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Capítulo 5: Teoria do Campo Autoconsistente (SCF) 


5.1 Introdução 


A teoria do campo autoconsistente (self-consistent field, SCF) para sistemas 
atômicos e moleculares, na formulação de Slater (1), se baseia no modelo seguinte: em vez 
de considerar todas as interações instantâneas entre os pares de elétrons do sistema, supõe- 
se que o potencial atuando sobre cada elétron é devido 1) aos núcleos e 2) à distribuição de 
carga média dos outros elétrons. 

Para os átomos com camadas fechadas esta distribuição é esférica, e o modelo é de 
campo central (1/r), logo um campo atuando sobre um elétron em ғ;, depende então da 


distribuição de carga média dos outros elétrons, que tem a seguinte forma: 
T ж 1 
Vrepulsão (4) = f Y (12,.,2N) s Y... 2N)d T..dTw [1] 
y 


ou seja, que está definido em termos da autofunção que é justamente a variável. O 
problema se resolve utilizando um método iterativo: postula-se uma certa distribuição 
eletrônica que permite calcular o potencial aproximado, constroe-se então o Hamiltoniano 
aproximado е acham-se as suas soluções. A nova V é supostamente melhor que a primeira; 
o campo construído a partir dela dá origem a um novo H cujas soluções dão origem à um 
novo campo, е o processo se repete até que o campo seja autoconsistente, ou seja, até que о 
potencial calculado em duas interações consecutivas seja constante, dentro da margem de 
precisão pré estabelecida. 

A teoria SCF, nesta formulação conduz às equações de Hartree-Fock-Slater, que são 
a base de métodos usadas especialmente em física do estado sólido, e popularizadas em 
química. 

Existem porém várias teorias SCF: a que os químicos usam geralmente é a de 
Hatree-Fock-Roothaan (HFR) (2), que é uma forma aproximada da teoria de Hartree-Fock 
(HF). Ambas as teorias se baseiam num modelo puramente matemático: para o caso de 
sistemas de camadas fechadas, consiste em representar o estado do sistema por uma função 
Y escrita como um único determinante de Slater em termos de um conjunto de funções de 
base monoeletrônicas Ivi J: sendo que as funções Y; são otimizadas utilizando o princípio 


variacional. No modelo de НЕ, as Ҹ;ѕдо variadas sem restrições, as equações são 


. . Ж x -— HF 
resolvidas numericamente, e as solugóes sáo os orbitais de Hartree-Fock, 52) E a 


energia calculada formado um determinante de Slater com estes orbitais, é a energia 
mínima que é possível conseguir com uma função de onda deste tipo. 


Na prática, as variações efetuadas em (V; } estão em vários sentidos, e os orbitais 
ш) não são estritamente calculados. No modelo de HFR, as equações não são 


resolvidas numericamente, e os orbitais são expandidos como combinações lineares de 


funções de certos conjuntos finitos, v m que podem ser fungóes gaussianas, orbitais 


atómicas, etc.: 


m 
(HFR) . (i) 
ч. =o Ф, [2] 
и 
a y I o (HER) x ada 
onde m é a dimensáo do conjunto de base d u ). Os orbitais Ф; são aproximações aos 
y (HP). 
j : 
"A (HFR) (HF) 
limite" =P; [3] 
т-»о 
No caso de moléculas, os orbitais moleculares (МО) são expandidos como combinação 
linear de orbitais atômicos (LCAO), e o modelo chama-se de SCF-LCAO-MO. A 
interpretação física das equações resultantes do modelo matemático acima não é tão 
simples quanto a do modelo de Slater: além da interação de Coulomb entre os elétrons, o 


potencial contém termos, chamados de troca ou de permuta (exchange) que surgem devido 
à antissimetria da função de onda. 


5.2 Valor médio da Energia 


Dada uma função de onda Total eletrônica para um sistema de 2N elétrons е М 
núcleos, escrita como determinante de Slater em termos de um conjunto de 2N spin- 
orbitais, 


Total í 
+ =A есіп! 5 Y Spin | (41 
Onde А é o operador antissimetrizador. А expressão para о valor médio da energia é: 


т Ге mar Р 
)= [vivaz P. 


Sabendo que a função de base é antissimétrica e o operador Hamiltoniano 


2N i> M 24 2N2N | 
Ер 6 
112 ATA ixj 


lembrando que este é um Hamiltoniano especial (a interação de um elétron com o restante 
do сагосо), a Primeira parte entre colchete está constituída por operador de um elétron, os 


termos de energia cinética dos elétrons е a porção da energia potencial relativa à interação 
dos elétrons com os núcleos: 


2N 
Hi => h(i) [7] 
onde | 
MES у, уй [8] 
i 2 ! A ҮД; 


Freqüentemente se utilizam aproximações nas quais os elétrons das camadas 
internas não são considerados individualmente; o efeito dos elétrons internos é introduzido 
no Hamiltoniano A(1) o qual contém a energia potencial do caroço molecular: 


2N 
Н} = x aa (D [9] 
1 
а segunda parte do Hamiltoniano: 
2N2N 1 
н УУ По 
i< j U 


e depende de pares de elétrons. 
Para calcular <Е> é conveniente escrever o operador antissimetrizador А em termos 


do operador permutação, P (umas vez que as funções de base são antissimétricas): 


me | Ps 
am? Dp [11] 


O fator (-1)” vale +1 segundo que a permutação seja par ou ímpar, e a somatória é sobre 
todas as possíveis permutações de 2N elementos, ou seja, que tem 2N! termos. O operador 


Ê permuta as coordenadas dos elétrons na função sobre a qual atua, assim, por exemplo: 
Рур MANO) Ge) Fo (4) 8(4)}= 


[12] 
= (DDD OLOY DF DLD) 
para determinar se a permutação é par ou ímpar, compara-se os subíndices, neste caso 342, 
com a ordem crescente normal, 234, e obtém-se o número de transposições necessárias para 
transformar uma seqúéncia na outra. No caso considerado, transpondo o número dois em 
duas casas para a esquerda obtemos a seqüéncia 234; como são necessárias duas 
transposições, a permutação é par, e (-1)' = (-1) = +1. 
Como a função Y [4] está normalizada, o valor médio <E> pode ser escrito como: 


(E)= e pP" f Piet Qe v oye Nm + Ho] 
[13] 


Б lv, ЕСЕ ый жы 


5.3 O Valor médio de H, 


Consideramos primeiro o integral sobre Hi, que chamamos de (H1). Como os 


elétrons são indistinguíveis, e como Y é antissimétrica, as contribuições de todos os /(1) 
são iguais; assim: 


2N 
= | ub ho var =2N[ Y n(iyYdr [14] 


Substituindo Y em [14]: 


НИ == aa DPP [Р OOR IO yaQ. ума) 
| 
h(1)P' lv, (Da (D, DADY, (3)aQ).. V y (20) 8QN)lr dradr3..dr) N 

15] 


(2 М") 2М!-1) 


а expressão [15] é uma soma de integrais múltiplas. E preciso analisar 


separadamente dois casos: 
(DP z P 

Suponhamos por exemplo que P é a identidade, e Р- В, , OU seja que Р permuta 
os elétrons 1 e 2. A integral múltipla correspondente é: 


ІҢ vr (e (DT Q)4Q).. y QN)BQN)RO)Y, (DADY (MIO. Y y (QN) BON dr dro d ray 


= Y DADO MEM Y IAS, Dadra [...[ Fx (LN) BON) y QN)BON)droN 
[16] 


O operador h(1) só depende das coordenadas espaciais do elétron 1, e não das suas 
coordenadas de spin; assim, a primeira integral do membro da direita da equação [16] pode 
ser fatorizada ainda em: 


ON aav, | 2 DAMA, 


е, como as funções de spin são ortogonais, а segunda integral é zero. O mesmo acontece 
com a integral sobre elétron 2. quanto às outras 2N-2 integrais, são todas de normalização, 
ou seja que são iguais a 1. 

Um caso levemente diferente ocorre se permutamos os elétrons 1 e 3 por exemplo. 
Temos: 
Г. Гэт DEDE DLDR Ga). y CN) AN) 
Y, (3)а(3)Ҹ DADY (а)... y QNBQN dr dro drs...d roy 
= [ет DADY, aa оаа (жү 050), DLAT |2 (3)% Gydvs 


ГезаолаҺ |... Y NN BON? y QN)BO.N)dro y 
Observamos que mesmo se as integrais de spin não são zero, mas a integral 


ж 
[92 G)av; 
se anula pois Y, e Y, são ortogonais. 
Assim, se P z P', todas as integrais sobre o operador Л(1) (ou, em geral, sobre 
qualquer operador de um elétron) são iguais a zero. 


(2)Р = Р' 


Obtém-se 2М! integrais com Pe Р iguais; para todas elas, o fator (pr é 


(=D 2P yle diferença reside no spin-orbital associado ao elétron 1. são (2N-1)! Termos 
para os quais 1 está associado a um certo spin-orbital V;o ou V; : 


[wi OY, 04 [Ге олов олвол со 
= [wi ORO; (а 


pois é possível fazer (23-1)! Permutações dos 23-1 elétrons, deixando o elétron 1 


associado ao spin orbital V; %- 


1 a 
tn) ON DID? [vi ORO; (Dav, 


N 


i 


о fator 2 aparece devido a que o mesmo orbital V; está associado às duas funções de spin 


a e В. 
5.4 O Valor médio de Н» 


O valor médio de (H 2) pode ser calculado de maneira análoga. Neste caso, para a 


função Y antissimétrica, a dupla soma: 


2» 
ie j 
pode ser substituída por 
2N(2N -1) 1 
2 n2 


_ 1 Р+Р' pap" T 
(нә)- збу 3j 1,2, CD f af 27 (Da)... y ONSON) 
1 


"12 


[18] 
P b. (Da)... CN)BCN)nidra day 


Consideremos separadamente os trés casos seguintes P= P', P = P' em mais de 
» 


dois elétrons, e P z Р' em dois elétrons somente. 
()Р= Ё 


Quando Ê= P', podemos fazer (23-2)! Permutações nas quais os elétrons 1 e 2 
estão associados aos mesmos spin-orbitais Fls e Tj %- Se os dois orbitais Y; e V; 


são diferentes (iz j), estes podem estar associados a quatro possíveis combinações das 
funções de spin: а(Йо(2); 5(1)8(2); a22); Pa), e para cada integral 


* * 1 
Jj = |: (09; Q — DF, Ojdr dr, [19] 
Ha 
На 4 (2N-2)! Termos. O integral J;; é chamada de integral de Coulomb. Portanto, se i=j, 


ou seja, se as funções espaciais do dois elétrons são as mesmas, as funções de spin devem 
ser diferentes, sendo que as duas possibilidades são: а()/ (2); Pal). 


Ju = [wi wi 0) vw, Qv, (244; [20] 
"n2 


(2) P = P' em mais de dois elétrons 


Se P difere de P' em mais de dois elétrons, é fácil demonstrar que as integrais valem zero. 


(3) P  P' em dois elétrons 

Se Р difere de P' em só dois elétrons, porém as funções têm apenas dois spin-orbitais 
invertidos (ou seja duas fileiras no determinante de Slater de uma das y). Novamente é 
preciso considerar vários casos: 

(i) se ambos spin-orbitais tem a mesma parte espacial, V; =P ;, a integral é do tipo 


[vi (Dav; 240) PA DEDE; Q)e (Dd ride» 


que é igual a zero porque as fungóes de spin sáo ortogonais. 


(ii) se as partes espaciais são diferentes, e os spin também, 
* " 1 
Ге (Цейуғ; А O E у DEDF;C)aldrdr, 
12 
a qual também é igual а zero por causa da ortogonalidade das funções de spin. 


(iii) se as partes espaciais são diferentes mas as de spin são iguais, a integral é, por exemplo 
* * 1 
[+ aY; (2)0(2) E Y,0Ma()Y,()a(2)d7,d7) [21] 
12 
que se reduz а 


K; = [BOPE О) тато [22] 
AV) 


a integral К; se denomina integral de troca (exchange). Se aos orbitais Ф; e V; da 


equação [22] se associam funções de spin В, o resultado será o mesmo. Ainda para cada 
para de orbitais (Ҹ;,Ҹ;) há 2(QN-2)! permutações dos outros spin-orbitais. Assim, a 


contribuição de cada par (*P;,W;) à energia total devido às integrais de troca é de 


2(2N — 2)K;. Finalmente, voltando à equação [18] temos: 


N NN 


i үжі 


Os Kj; aparecem negativos porque correspondem à uma permutação impar. 


5.5 Cálculo da Energia Média 


Juntando os dois valores médios [17] e [23] temos: 


N N NN 
(E)-29 hi9 J; +) (J; - Kj) [24] 
i i i jei 
esta expressáo pode ser simplificada observando que 
* * 1 
Ji =K; = Із; EA) A 0 206 1878 
12 


de modo que 


N NN 
(E)=25 + УУ (023; Ky) [25] 
i i j 


Cada termo da primeira somatória na equação [25], 24%, corresponde à integração entre um 
elétron no orbital V; e o núcleo; o fator 2 é devido a que dois elétrons ocupam cada orbital 


'P,. Os termos 
N 
озу -ку) 
j 


representam a repulsão entre um elétron no orbital V; e os outros elétrons. E interessante 


definir um conjunto de energias orbitais E; para um elétron ocupando о orbital \Р; 


N 
Е, = +$ (J; - Kj) [26] 
j 
Se fizermos a soma destas orbitais sobre todos os elétrons do sistema, teremos: 


N N NN 
25 E; =2У hy +25 > QJ - Kj) [27] 
i i i j 
comparando esta expressão со a da energia total «E^ da equação [25], observamos que 


N NN 
(E)=25 E; +) QJ j - Kj) [28] 
i i j 


o qual é facilmente justificado já que na energia total, a repulsáo entre cada par de elétrons 
em Ф; e Y; deve ser contada só uma vez. 


5.6 As Equações de Hatree-Fock. 


Na seção anterior vimos que, dada uma função de onda aproximada, escrita como 
um determinante de Slater em termos de um conjunto de funções de um elétron {Ф;}, а 
expressão geral para a energia do sistema é: 


N NN 
E=29 h; 429 3. QJ - Kj) 
i i j 


podemos aplicar o método variacional para achar as melhores funções de base V;, isto é, 
aquelas tais que se formamos V como um único determinante de Slater na base das V ;, а 


energia é mínima: 
óE = 0 


А variação do conjunto teil porém está sujeita à restrição de que continue sendo um 
conjunto ortonormal 


Sy = 0 se іж/ 
] se i=j 
ou seja, que: 
óS; -0 


O método que permite variar E mantendo a ortonormalidade da base é um método 
matemático bem conhecido, o dos multiplicadores de Lagrande. Formamos uma função G: 


G=(E)- 22 455 


onde os 4; são os multiplicadores de Lagrange, por enquanto, parámetros indeterminados, 


e procuramos o mínimo na função С: 


óG -óE - > > A568; =0 
i j 
Substituindo E por seu valor (equação [25]) 


NN 
25 hj; +) 2 061;-дк,)-У2,4;05; =0 [29] 
i i j i j 


Os termos da primeira soma na equação anterior: 
Shii = o v; DADP; Mar, = | (#Р; DADY, (dri + [Y (010% (Mar, 


podem ser reescritos como: 


Ой; = СА (DA(O)V; (Dd v, + complexo conjugado 
Analogamente: 
* * 1 * * 1 
à = f OP (УФ ;Q) — V, (D, Odrdr + f Vi MO), DP) Odridr + сс. 
12 12 


ou utilizando o operador integral de Coulomb: 


R A 1 
3,0) 2E os (2)4т» [30] 


Ay; = | (9, (1))3 PE, (Ddr + [ (FME; (Ddr + cc. 
Para as integrais de troca 


Ж = [C7 OR OY, Ddr + [OK Qj (Ddr + ce. 


e definimos o operador integral de troca, K j (D. como um operador que tem a seguinte 


propriedade: 


К ,(Dv,()- | pe o (2)аг› | ¡0 [31] 


Finalmente, a variagáo das integrais de recobrimento é dada por 
ж 
Ty = [O7 MF (Dari + сс. 


Juntando todos os termos e pondo em evidência as diferenciais: 


N N 
Ж f (SP (14 2h Y; (1) + 25,05 ¡MK 00) Aj Y; Dc tcc. -0 
i J J 
O fator 2 que aparece na frente da primeira soma sobre j é devido à que tanto Jij como Jj 
contribuem com um termo idêntico,e analogamente К; e Kj. Como as variações das 
funções e dos seus complexos conjugados são independentes, cada termo da soma em i 
deve ser identicamente igual a zero. Assim: 


N N 
20)+2 Q3 (0) -K;() V; - 345 V;() 20, i=12,.N 
j j 


substituindo, por conveniéncia, 


e dividindo tudo por dois, temos: 
N _ А М 
ADS QJ, 0) -K A W,0)=> E;V;(0)-0, і-і2,.м [32] 
j j 


A Expressão contida entre os colchetes se denomina operador de Fock, Ë : 


A N ^ ^ 
Е -h() + $,QJ (D-K ;()) [33] 
j 


temos então um conjunto de equações diferenciais de um elétron: 


N 
Ғғ ()- ОХ Ej PD, 151,2,...N [34] 
j 
Vemos imediatamente que se todos os Ej; com i= j fossem iguais a zero, a equação [34] 


seria uma equação de autovalores (devemos diagonalizar Ej para adquirirmos Ер). Isto 


pode ser conseguido mudando-se o conjunto de base hi }. Com efeito, os orbitais Y; não 


são únicos. Qualquer combinação linear, 
j 


também é autofunção do operador É, se os coeficientes Tj (transformação unitária) são 
tais que: 
ж 
k 


uma transformação que obedece o requisito [35] é uma transformação unitária. Em forma 
diferencial é 

T'T=1 [36] 
onde Т” é a matriz adjunta de T, ou seja, a matriz transposta conjugada. A equação [34] 
pode ser escrita em forma matricial: 

FY = EY [37] 


Multipliquemos ambos membros desta equação рог T pela esquerda, e insiramos a matriz 


1, ou, o que é equivalente, o produto TYT: 


ТЁТ*ТФ = TET* T 
Agrupando os fatores, temos: 
(TÉT?)TY = (TET* )TY 
ой 
ЭРЕ 


que é análoga à equação [37]. Em particular podemos escolher uma matriz T tal que E'seja 
diagonal. Neste caso, teremos uma série de equações de autovalores. 


PY = ENE, i=1,2,...N [38] 


para o operador Ё". É possível demonstrar que o operador É" na base (F';! é exatamente 


análogo ao operador É na base teil de modo que as equagóes [38] representam o mesmo 


problema que as equações [34]. As equações [38] são as equações de Hartree-Fock, e os 
orbitais {\Р'; } são os orbitais canónicos de Hartree-Fock. 


Podemos escrever a equação [38] diretamente 


Ғе; - EY; 1=1,2,...N 
com 
F=kH0D+03;-K)) 
j 
Notamos o Valor E;: 


E; = [viFw;dr = Е; 


| %; h() + QJ; -K y) Vidc =hj +X QJ; - Kj) 
J J 
é exatamente igual a energia total E; definida pela equação [26]. Estas são chamadas de 


potencial de Koopmans ou potenciais verticais. 
A resolução das equações de Hartree-Fock se faz da maneira seguinte. Como o 


operador de Fock, É , contém os operadores integrais J je K j definidos em termos de um 
conjunto de orbitais de base teil inicia-se o processo postulando uma forma para o 


conjunto hei por exemplo, os y O podem ser autofungóes do operador h(1) 


hw (т) = gy (0) 


com essas funções calcula-se o operador de Fock е resolvem-se as equações de Hartree- 
Fock. As soluções е) formam um novo conjunto de base по qual pode se definir 
novamente o operador de Fock. O processo continua até que os orbitais permaneçam 
invariantes sob novas iterações, dentro da margem de precisão requerida. Se os yO (1) 


estão bem escolhidos geralmente o processo converge. 

É interessante discutir o significado físico do operador de Fock: ele representa um 
Hamiltoniano efetivo de um elétron tal que a energia potencial do campo no qual o elétron 
se encontra é constituída das seguintes partes: 

(1) A energia de interação com os núcleos 

(2) A energia de repulsão com todos os elétrons de spin oposto ao do elétron 

considerado. 2: J; 

J 

(3) А energia de interação com todos os elétrons de igual spin ао do elétron 

considerado: 


Db; -&,) 


que é a interação (2) devido a que, para uma função antissimétrica, dois elétrons com o 
mesmo spin não podem ocupar o mesmo orbital e consequentemente nunca estão muito 
perto um do outro; assim a sua energia de repulsão é menor 
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Capítulo 6: Método de Hartree-Fock-Roothaan. 


6.1 Introdução 


Para sistemas com muitos elétrons a resolução das equações de Hatree-Fock pode 
ser feita de maneira aproximada pelo método de Hatree-Fock-Roothaan no qual os orbitais 


de Hartree-Fock fr}, são aproximados como combinações lineares de funções de base 
iz x: 


>= Усх, [1] 
H“ 


para moléculas, os \Р; são orbitais moleculares e os y, são orbitais atómicos. O método 
(i) 


variacional é aplicado aos parâmetros lineares с. 


Calculemos primeiro as integrais Л, J; e К da expressão abaixo, em termos das 


funções Ku] 
N N NN 
2» E; =2) h; 429 > QJ y - Kj) [2] 
i i i j 


sendo assim obtemos, 


hi = [97 020% = Y eO e | Олут =D cha — [3] 
uv uv 


* * 1 De (y (DE (GF ж ж 1 
1; = LOGO, DEDO -X X XX DO | 0,0) — xis dnde = 
ho и у А © "2 
= YU e e S) (uA | vo) 

ШАС 


[4] 


: 1 
onde o símbolo (uA|vo) se refere à integral sobre — tal que as funções da esquerda estão 
"12 


associadas ao elétron 1 e as da direita, ao elétron 2. Analogamente: 


Kj = y u t (272 0 сід (uo | vA) [5] 
vio 


Substituindo na equação abaixo 


N N NN 
i i i j 
temos: 


N m NN А 
E= Y cn УУ уусу e ee Aa vo) (ua vay] 16) 
i i j uÃo 


onde é conveniente definir: 
N 
E (i)* (1) 
Puy = 22, Cu С 
1 


Substituindo па expressão de E: 


E-Y Ры, +5 DP sa Pro пано) - ue A)| 
ШУ 


2 


) 


procurando agora os melhores valores das constantes D variando a energia com respeito 


a essas constantes. O procedimento e inteiramente análogo ao da secção anterior, e o 


resultado e, 


"TEE D N 
Oe os +A» cDDP (ua |vo)- (uc | v2)] = 2.692 649 
Y j j у 


onde y =1,2,... 


fazendo uma transforma-c”ao unitária da base para diagonalizar a matriz e, e definido 
x 


Fw = + Y Pro QA | vo) - (uo | v4)] 


vo 


obtemos 


Y us – ES Je? =0 [7] 


que são as equações de Roothaan. Em notação matricial, torna-se 
FC = SCE [8] 


A grande vantagem pratica destas equações e que elas são simplesmente equações 


algébricas e não mais diferencias. 


Capítulo 7: Teoria da Perturbação. 


7.1 Introdução 


Agora nós iremos discutir um outro método de aproximação em Química Quântica, 
a teoria da Perturbação ou Teoria da Perturbação de  Rayleigh-Schródinger. 
Consideramos um sistema o qual a equação de Schródinger é 


Н|\Р; у= Е;|\Р;) [1] 
e para um estado estacionário 
onde H 9) é Hamiltoniano de um sistema não perturbado, ou seja um sistema na ausência 
de radiação eletromagnética ou correlação interna. Quando há presença de uma perturbação 


o Hamiltoniano que descreve a perturbação é representado por H”. O Hamiltoniano 
Perturbado H e um não perturbado HO é 


H=H-HO [3] 
H=H04H' [4] 


onde HO > gr. Quando dizemos que H (0 5 H ', queremos dizer que (н 2i » (H ” Я 


isto é, estamos comparando os valores esperados dos dois operadores H Оен. 

O objetivo é relacionar as autofunções е autovalores do sistema perturbado às 
autofunções e autovalores conhecidos do sistema não perturbado. Para isto, imaginamos 
que a perturbação é aplicada de maneira a variar continuamente o sistema do estado não 
perturbado ao estado perturbado. Matematicamente, isto corresponde à introdução de um 
parámetro A по Hamiltoniano, modo que 


H-HO +15 [5] 


com À variando entre O e 1. Quando 4 =0, o sistema é não perturbado e quando 4-1 o 
sistema é perturbado. Assim temos: 


HO + Amy; = E; v, [6] 


desde que o Hamiltoniano depende do parámetro A, tanto a autofunção V ;, como o 
autovalor E;, devem depender também de 4 e podemos expandi-lo em séries de Taylor de 


potencias de 4. 


y; PO D OL +...+л*ч® [7] 
Е, = EO + Ag + ED +. + АЕ [8] 


Onde k é a ordem da função e energia. Substituindo [7] e [8] em [6], temos: 


(HO + дне + л? + 290) + 0 = 


[9] 
(EO + Ag(D + REPOSO ++ Me) 


HOPO + А(Н'Ф{? - HOW(O) 4 22 (H OO) - pO), = 


[10] 
ЕО + А(Е GO + ENG, +4° (gw + EO + ЕО +... 
Igualando as famílias de mesma potência 
HOYO = pO» [11] 
np (0 0 1 1 0 0 1 
AER + gi p 0) = (gp + (D) 
np (Ü 0 1 1 0 0 1 
HYO 4 HOYO EDO + p -0 [12] 
22 ( HOYO Т av, = 22 ( Ep 4 Е0ҹ%ҹ(0 + ЕО) 
н Ow 02) к Hy (0 = EO (0) 4 E Ow (l) 4 Ey О) -0 [13] 


Agora devemos agrupar as expressões [11], [12] e [13] 

(HOPO -gO O) Дн-0 yp O 4 (gi _ gy (D |+ 22(.)+..=0 [15] 

se V; e Е; são funções contínuas de A, рага que a equação [15] seja verdadeira para 
qualquer valor de 4, os coeficientes de cada potência de 4 devem ser iguais a zero. As 
equações [11] e [12] são equações de ordem zero e ordem um respectivamente. E como as 


soluções de ordem zero são conhecidas elas podem ser utilizadas se obter as equações de 
ordem um. 


7.2 Correção de ordem um para a energia 


Para obter a energia reescrevemos a equação [12] 


н'Ф®) HO E) -gy0 
1 1 1 1 1 
multiplicando por y OP ambos os termos e integrando 
ғ” HYO E ә” (Н (0) _ go yw = y OP EO (0) 


[um | HO -gO | №0) 2 ЕШ | н 


v - gi 


O primeiro termo é zero porque H é um operador Hermitiano. Então, a correção de ordem 
um para a energia do nível i é: 


Е)-|ре УнФФОлтену [16] 
ou 

ED = © | HO | v) -H 

ED = nm | HO | yO ) т us | HO | y ) 

: [17] 


E => БІР PS өз, 


m=1 
7.3 Correção de ordem um para a função de onda 


Para obtermos a correção de ordem um para uma função de onda, expandimos yO 


em um conjunto ortonormal completo de funções de ordem zero, peo E 


Уа бийк Re |ui [18] 
j 


Substituindo [18] em [16] temos: 


Y a, 10 _ 609900) =(E0 -NPO 
j 
ou 


Na (E EN ES - my [19] 
J 


multiplicando а equação [19] por yo e integrando, temos: 


La (EO ED |p) = EP (O O) (e A) ро 

j 
Efetuando a somatória, somente os termos com k = | permanecem, devido à 
ortonormalidade: 


AE EE RO |в) (нере 


m [21] 
temos agora dois casos a considerar: para К-і, obtém-se a equação [16]. 
КОШЕСІ 
— — — 


ЕФ = (w | н 


um Ж DEP um m [16] 
Рага k +i, о primeiro termo à direita da igualdade é zero e a equação [21] torna-se: 


a (E -Е@®у=-(Ф®|н'|®®) 
arranjando 


= RS cpm [22] 


a equação [22] determina todos os coeficientes а, exceto а,, o coeficiente de YO. 


ii? 


Aplicando a condição de ortonormalização, se determina а; = 0. Os resultados finais para 


E, e Y, corretos até a primeira ordem são os seguintes: 


ED = E) 44H , [23] 


(Erro о pa 


YO au 0 41) 
г d (0) (0) Í 
EO — gt 


kzi 


7.4 Tratamento de um sistema de spins pela Teoria da Perturbacáo 


Quando estudamos um ensemble de spins através de perturbações, realizamos um 
deslocamento de equilíbrio do sistema de spins e podemos fazer transições entre níveis de 
energias que ocasionam mudanças na magnetização na magnetização líquida. Esta por sua 
vez nos dá o sianl de RMN. 


Considerando um sistema de dois níveis representados por |а) е |b), com energias 
Еа e Eb, sujeito a uma perturbação V(t)-Vf(t), onde V é um operador independente do 
tempo e f(t) é um fator numérico flutuante. O operador V, tem com efeito a modulação nos 
níveis de energia Ea e Eb, o qual induz transições entre os níveis. 


V|E,) = Va FO) Es) 
V|E,) = Via (O Es) 


Temos para esse sistema o Hamiltoniano 
үт | Е, L 
VJ Q) E, 


е a função de onda que é solução da equação de Schródinger: 


м) = c oep- Ja) e exp E 
Е 


2 
ih— 


E,t 
a c, (t) exp( 19) 


Resolvendo o sistema abaixo: 


E, V,f@ een] 
Tz D oel 


с> E 
inae, (t) = Vie OSO — 


temos: 


E 
е) = V,,e, (t) f (t) expli EED 


Supondo que no instante t=0 o sistema esteja no equilíbrio, ou seja, no seu estado ja) , em 


que ca(0)=1 e сь(0)=0. A integração da equação acima fornece o valor da correção em сь 
num tempo / posterior: 


"E [^ ^: f(r)ex wi et - 22 


ot) - -ie | ЕНШІ Ga ын a 


A probabilidade de que ocorra uma transição do estado ja) para o estado |b) no tempo ѓ ё 
dada por: 
2 
P(a,b) = lc, (2) 


Um caso em que P(a,b) cresce linearmente com o tempo é aquele em que ҚО) é uma 
função periódica do tipo: 


f (t) = 2 cos(at) 


-Е 
c, (t) = ғы, cos(at) ew; (Es : al Et, 
(ЕЬ- Ea) Bus | do 
escrevendo @„ = s (ou seja, hou, é a diferença de energia entre os níveis) 
Cy (t) = | 2 cos(at) expli о! 


е вепдо: 
exp(iat) + exp(—ioxt) 


cos(ot) = > 


1 


c(t) = –і e f (exp(i ot) + exp(-i ot)) exp(ioy,t)dt 
0 


ү t 
c, (t) = —i 3 | (ехр@ ot) + exp(—i ot)) exp(i 0.0)! 
0 


er) - -i al | (explico + op r] + explitoy, on 


0 


Opa FO Opa = O 


a 


ce Kia ЕБЕ +o)t]- l. expli(o,, = @)t] zl 


Рага œ próximo de оъ podemos desprezar o primeiro termo da equação acima com 
relação ao segundo termo.A probabilidade P(a, b) = |с b (г? fica então igual a: 


P(a,b) = c, (T| «Pal estt, 0101, ton -01-1 


Oba = Ф Op, — Ф 


2 | | 
Лу. | ехо, — а) -оу) +1 
n (oy, о) 


P(a,b) = Voal | 2 — 2cos(@,, — t 
| a (Opa - 0) 
P(a,b) = LAE 1-сов(о,, - o)t 
uam "m 
| 2 Ба 
у Е зал (ды 24 
P(a,b) =— 
1 
(05, -о) 
sa © - y | 
2x V, 
P(a,b) - яуы) 2 


n MCN 
2 


P(a, 


Pm 
1 1 1 1 1 1 1 L 1 1 1 1 1 1 
Oba [0] 

probabilidade de transição em função da frequencia da perturbação (fig.3) 


Na figura 2 temos uma simulação da função probabilidade de transição P(a,b) entre dois níveis 
com energias Ea e Eb. Aproximando essa função por uma função delta de Dirac: 


Lxx, 
СЕ 0;х Z x 
? 0 
2z|V,, А 
LU тин - ot 


A primeira condição a ser satisfeita pela perturbação é que sua atuação forneça ao sistema 
uma energia hw=hv, onde у= (Eb-Ea)/h, ou seja, a energia fornecida ao sistema deve ser 


igual à diferença de energia entre os níveis. 
O termo independente do tempo da perturbação deve ser tal que | Vpa | 40. 


ЖЕСІ = 0 
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